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We review an exact WKB resolution method for the stationary ID Schro- 
dinger equation with a general polynomial potential. This contribution covers 
already published material: we supply a commented summary here, stressing 
a few aspects which were less highlighted before. As some of our earlier 
papers needed later corrections, we also recapitulate these here in footnotes 
in the bibliography. The latter, not meant to be exhaustive, focuses on 
narrowly relevant works plus some broader ones that do offer more extensive 
bibliographies. 



1 General setting 

We will treat a Schrodinger equation on the real line with a polynomial poten- 
tial function V (real, normalized as V{q) = + Vi q^^^ + ■ ■ ■ + Vn-i q)'- 

(-^ + [n?) + A])^(g) = o, (1) 
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versite Paris 7, F-75251 Paris CEDEX 05, France. 



1 



(a general Sturm-Liouville problem). It is convenient to pose the problem on 
a half-line, initially q G [0, +00); exact treatments however rely on analytical 
continuations to the whole complex domain: g G C, [v, A) G (where 
V = {vi, . . .,vn~i))- 

1.1 Outline of results [,47j-[50| 

The traditional view about eq.(Q]) (in any dimension) is that 

- iV = 2 is exactly solvable: the harmonic oscillator, "the nice" case; 

- no other N is solvable (generically) : 

• = 1 gives the Airy equation (more transcendental than the higher- 
degree iV = 2, a paradox!) 

• > 3 yield the anharmonic oscillators, the ground of choice for all 
approximation methods (perturbative, semiclassical . . . ). 

We see eq.(^ instead as an exactly solvable problem for any A^, 
thanks to an exact WKB method with two key ingredients: 

- a semiclassical analysis exploiting zeta-regularization |^ I29j: 

- exact quantization conditions that are Bohr-Sommerfeld-like but selfcon- 
sistent, and akin to the Bethe Ansatz equations that solve many exact 
models in 2D statistical mechanics ||2ri EU] or quantum field theory |7| . 

Moreover, degrees N and 4/iV behave similarly in this framework (a qual- 
itative duality) |1B] : 

• = 1 is transcendental like A^ = 4, clarifying the above paradox; 

• N = 2 (selfdual) becomes completely explicitly solvable as a degen- 
erate, rather pathological, case. A countable family of similarly solvable 
cases is actually obtained [SUl Sec. 4]: all the generalized eigenvalue problems 




+ (g^ + Ag 2 1)] ■?/)(g) = (for even degrees A^ > 0). 



1.2 A few early sources 

Exact asymptotics (which uses drastically divergent representations to per- 
form fully exact calculations) is a paradoxical idea; its advent rests on var- 
ious seminal developments which mainly took place in the '70s (restricting 
to works we felt as directly influential; a comprehensive bibliography would 
need much more space): 

- Balian and Bloch's representation of quantum mechanics, where the lat- 
ter can be exactly rebuilt in principle just from the classical trajectories. 
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provided these are taken in their full complex extension jH]; 

- the exact geometrical content of high-frequency asymptotics: singularities, 
wave front sets, and their propagation properties ||261 [13 122] ; (reviewed in 
EH chap. VIII]) 

- the same in the analytic (not C°°) category (a key step to exactness): ram- 
ified Cauchy problem [S2], analytic pseudo-differential operators 11 , hyper- 
and micro-functions localized Fourier transforms ^21; (reviewed in 

- Borel resummations of asymptotic series [201; large-order perturbation the- 
ory P, tunneling (e.g., [23 [HI 133] ) , instantons and Zinn- Justin's conjectures 

- seemingly unrelated direct monodromy calculations resulting in exact func- 
tional relations for Stokes multipliers |38j . 

The Balian-Bloch idea led to a first concrete exact WKB algorithm for ID 
problems like eq.(^, based on (generalized) Borel resummations [1^ [^ [T7| 
nHlEll2Hl! and best formalized in Ecalle's framework of resurgent functions 
[23; the explicit exact WKB results are indeed special instances of bridge 
equations in resurgence algebras (or resurgence equations). Still, in this 
approach, one basic resurgence theorem is not yet fully proved |181 thm 
1.2.1 and its Comment], and a fairly explicit solution stage has only been 
reached for homogeneous polynomial potentials, rather tortuously [^-|46j. 
Otherwise, many developments have been achieved in this and neighboring 
directions: proofs of Bender- Wu [301 [IH] and Zinn- Justin conjectures [T71 ITH]. 
hyperasymptotics (reviewed in [2H1 Part II]), results on Painleve functions 
(reviewed in [28^ Part III]). 

This presentation is about a more direct exact WKB method |47j-[SUj 
(surveyed in ^51] ), based upon functional relations a la Sibuya plus cru- 
cial inputs of zeta-regularization both at the quantum and classical levels. 
At its core, it needs just one exact Wronskian identity plus a structure equa- 
tion for a pair of spectral determinants, to generate the exact solution as 
a fixed point of an explicit system of nonlinear equations ("exact quanti- 
zation conditions", but also Bethe-Ansatz equations in an integrable-model 
interpretation [21] )• The Wronskian identity precisely subsumes the resur- 
gence relations found in the Borel-transform approach, whereas the newer 
exact quantization conditions (which achieve the solving task) have no clear 
counterpart in the Borel plane as far as we know. This direct approach also 
makes it easier to treat the general case at once, and it is now fully proven 
for homogeneous potentials [1]. 
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2 Basic ingredients 

2.1 The conjugate equations 

Just as a polynomial equation is better handled by treating all its conjugate 
roots together, the original differential equation (0) is to be supplemented by 
its conjugates, defined using 

''^ e-''^ V{e-''^/' q), = e"^^^ A, with <^ ; (2) 

there are L distinct conjugate equations labeled by £ G Z/LZ, where 

+ 2 generically 
y + 1 for an even polynomial V{q). 



2.2 Semiclassical asymptotics 

• WKB formulae for solutions of the Schrodinger equation jHE]: a solu- 
tion of eq.(0) which is recessive (= decaying) for q —>■ +oo is specified by the 
WKB behavior 

Uq) - n,(g)-i/2 exp -Tu.iqW, ^,(g)'=^'(V(g) + A)^/^ (4) 

JQ 

where Ilx{q) is the classical forbidden-region momentum {Ilx{q) dq is the 
Agmon metric). Reexpansions in descending powers of g — *■ -|-oo yield 

[V(^)+A]-+V2^^^^(,.^^^)^.-A^. (a=f, f-1, f-2, ...),(5) 

cr 

Mq) - e^g-^/^-'^-(^^)exp|-5^/3.„i(t;)^| (6) 

where I3a-i{v) =^ /3a— 1(0; v, A) (independent of A for cr > when > 2). 

Another solution, recessive for q +e~"^/^cxD, arises through the first 
conjugate equation: 

^A(g)'=^^S,(e-/^g). (7) 

As the asymptotic form (jHl) and its analogs for '^\{q) and their g-derivatives 
have overlapping sectors of validity, this explicit exact Wronskian follows: 

^Mi^'xiq) - ^A(5)^A(g) = 6^+^^'^' 2ie'^/^e'^^-('^)/2 0). (8) 
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• Spectral asymptotics (with E = —A): the Schrodinger operator over 
the whole real line with the confining potential V^di?]) has a purely discrete 

spectrum {-EfcjfceN {Ek T +00). Due to parity-splitting, 8^ =^ {E2k} is 
the Neumann, and £~ '= {E2k+i} the Dirichlet, spectrum (the boundary 
conditions being at g = 0). As a corollary of the WKB formula (jH) on eigen- 
functions, the eigenvalues obey the semiclassical Bohr-Sommerfeld condition 

S{Ek) ~ A; + i for integer k — >• +00, (9) 

in terms of the classical action function S{E) =^ §{p2+v{q)=E}^^^ which be- 
haves as b^E'^ for E +00 with = ^ + (the order). 

Further expansion of eq.(jni) to all orders in k~^ followed by reexpansion 
in descending powers of E results in a complete asymptotic {E +00) 
eigenvalue formula, of the form jlT] 

h^Ek^ ~ + I for integer k +00 (a = /i, fi — j^, ^ — j^, . . .) (10) 

a 

(the ba are polynomial in the {vj}j<{^~a)N, and also parity-dependent but 
not above a = —3/2 nor if V is an even polynomial). 

2.3 Spectral functions [143] 

• Spectral zeta functions (parity-split for later use): 

Z^{s, A) = J2 i^k + A)"' (Re s > fi), (and Z{s) = Z+{s) + Z-{s)) 

1 even 
odd 

(11) 

extend to meromorphic functions of s G C, regular at s = 0, all due to 
eq.dini). 

• Spectral determinants: formally meant to be 

D^{X) = " n i^ + Ek) ", (12) 

, oven 
^ odd 

they get rigorously specified through 

log {X) —dsZ^{s, X)s=o {zeta-regularization) , (13) 
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and further evaluated by limit formulae which we call structure equations 
(they are like Hadamard product formulae, but with no undetermined 
constants whatsoever, a crucial property here): 

logD^iX) = lim I V log(Efc + A) + i log{EK + X) 



(k.KZl) , (14) 

(counterterms:) — | baiEx)" log Ek — ^ 

{a>0} 

As corollaries for D (and D^): 

- D is an entire function in (A, v), of order fi in A; 

- logD has a large-A expansion of a severely constrained ("standard" or 
"canonical" ) form, or generalized Stirling expansion [221 • 



logZ^(A) ~ aiA(logA 



+ ao log A 



(15) 



-oo 



(degree 1) (degree 0) (degrees /i^ i 

where obviously no term can occur with degree > /i (growth restriction; e.g., 
ai = as soon as > 2); but above all, 

any extraneous power terms fen A" (n G N), including additive constants 
(6oA°), are banned outright (canonical constraint). 

For N = 2, eqs. lfT^ . resp. (fT^ . essentially restore classic results: the 
analytical continuation by the Euler-Maclaurin formula, resp. the Stirling 
expansion, for —dsC{s, X)s=o = log \/27r/r(A) (Lerch formula), with ((s. A) = 
the Hurwitz zeta function. Those two founding formulae of all asymptotics 
involve zeta-regularization, as the Lerch formula shows; by inference, WKB 
theory (a branch of asymptotics) also ought to be intimately tied to zeta- 
regularization, but this aspect traditionally goes unnoticed. 



3 Fundamental exact identities 

3.1 Canonical normalization of recessive solution 

The solution ^Jxi^q) described by eqs.(jH), (jHI to be recessive for q +oo, still 
awaits normalization. It proves useful to fix it also at g = +oo, as 

r+oo 

~ nA(g)-^/' exp / Ux{q') dq' {q ^ +cx)). (16) 

Jg 
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However, /+°° nA(g') dg' ("the A gmon distance from q to +00") diverges: we 
then wish to define this "improper action integral" as 



{V{q') + \f'^''Aq' 



(17) 



(analytical continuation from {Re s > /x} to s = 0), but s = is a singular 
point in general: a simple pole of residue N~^P_i{v) (cf. eq.®). Now: 

• /9_i(s; v) = quite frequently (for any polynomial V{q) of odd degree N; 
also for any even polynomial of degree N = 4M; for any V{q) = q^ , N 2); 
then s = becomes a regular point and the prescription (fTTjl is fine; 

• when /9-i(s; v) 7^ 0, a stronger regularization is suggested by the quantum- 
classical correspondence. "Classical spectral functions" are naturally defin- 
able, by mimicking the quantum formulae (fTT|) . (fT^ . as 

zds,x)= [ ^{p' + vm + ^r (18) 

logZ^cKA) = -9,Zd(s,A),=o. (19) 
These formulae together imply a classical counterpart to the formal eq. ()12|l . 

7-r(nkl) + A)^/Mg"=logDei(A); (20) 
using parity symmetry, this further splits to give the pair [201 Sec. 1.2.2] 

D±(A) = n,(0)±i/2exp "/;°°n,(g)dg "; (21) 

logi5^(A) are also the divergent parts of the large- A expansions (fl^j) for 
log L'^(A), so they are canonical as well. The idea is then to define °° ^x{q) dg 
to be I log Del (A) (and J^°° = Jq°° — Jq)- This in turn fixes the recessive 
solution, with a normalization e'' ^ 1 in eq.® (an anomaly factor, which 
would turn up elsewhere if not here): 

C ^ ±[-(2log2)P^,is;v) + ds{^^)]s=o • (22) 

As examples of the resulting improper action integrals: 

/;^(g^ + A)V2dg = -|7r-V2r(l + ^)r(-/i)A'^ {N ^ 2; /i = i + ^) 
j^°°{q^ + A)^''^ dg = —\ A(log A — 1) (anomalous when A 7^ 0) 

jt^iq^ + vq^l^Aq = -\v^'^. 

(23) 
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Remarks: 

- a "classical determinant" Dd (or D^) is not an entire function of A like D: 
the discrete spectrum of zeros of D classically becomes a continuous branch 
cut for Del- 

- our treatment of the normalization for ^ was initially flawed 

we erroneously prejudged the natural answer to remain C = 0, so we had tem- 
porary inconsistencies of normalization in the anomalous case; these affect 
eqs. (15-19), (26) in |17j and Sees. 1.1, 2.1, 2.2 in [SUl, but cancel themselves 
spontaneously thereafter; the error is repaired from jlH] onwards. 

3.2 Basic exact identities 

Under the canonical normalization for the recessive solution, and with ' = 

D+(A) = -^1(0), resp. D-{X) = MO)- (24) 

The idea, assuming > 2, is first to prove the logarithmic A-derivatives of 
eqs.(j21I) by analytically computing traces of the resolvent kernel, and to A- 
integrate back: then, no integration constants can reappear precisely because 
all four logarithms have the canonical large-A behavior ()15p. in particular 

f,iO)--Dt^W. resp. ^,(0) ~ Z},l(A) (25) 

(from eqs. (PT|) and (fT^ at g = 0, also valid for A +oo). In depth, eqs.(j21I) 
arose by controlling all integration constants without exception, from eq. ()14j] 
onwards. 

The basic exact identities can now be substituted into the exphcitly 
known Wronskian (jS)) computed at g = 0, resulting in a functional relation 
between spectral determinants or "Wronskian identity" , also fundamental: 

e+^^/^DW+(e-'^A)D-(A)-e-^^/^D+(A)DW-(e-'^A) = 2 ie+^^^-i^^^)/^ _ (26) 

4 Exact quantization conditions 

At first glance, the exact functional relations (j^Bj) look grossly underdeter- 
mined; like the original Wronskian formula (jHl), they amount to one equation 
for two unknown functions. The exact information here seems to present a 
"missing link" throughout. This is however a delusion, because those un- 
known functions have a truly special form here, embodied in the structure 
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equations (fT^ (themselves consequences of zeta-regularization), and this pre- 
cisely fills that information gap. 



4.1 Degenerate cases 

• = 2 : [— ^ + {q'^ + A)] ij{q) = (harmonic oscillator). 
This is the most singular case: = A/2 ^ 0, and it has A-dependence. 
Now, (fi = IT makes eq. ()26p degenerate: for A real, it separates into real and 
imaginary parts, which recombine as D^{X)D~{—X) = 2 cos |(A — 1). The 
zeros are then obvious, and they clearly split between -D^(A) or D^(— A) 
according to their sign; thereupon, the structure equations (fT^ yield the 
explicit solutions 



D^(X) - ^''^'^^ D-(X) - '''""^ (27) 
D (A) - D (A) - (27) 



• even: the generalized eigenvalue problems [— ^ + (g^+A g 2^ ^)] ^ip^q) = 
are also exactly solvable (SUl Sec. 4]: setting z/ = eqs.(OU|) and (fT^ like- 
wise lead to the relevant Neumann and Dirichlet determinants for q G [0, +C)o), 
as 

2-A/^(4z/)^(^+^+V^r(-2z/) T^-(^\- 2-^/^(4z/)''(^-i)+i/2 r(2z/) 



r(i^(A - 1) + 1/2) ' ' r(KA + i) + i/2) 

(28) 



whereas over the whole real g-axis. 



dptf-di + r«^ + A«f-i^l-/^^(^)^^(^) ^^^ = 2 (mod 4) 

(let [ + +A g 2 jj _ | ^^.^ if = (mod 4). 



(29) 

Remark: actually, the potentials (g^ + Aga i) are supersymmetric at the 
generalized eigenvalues A (see |50j and the footnote under that reference). 

4.2 The generic situation 

The preceding treatment does not generalize straightforwardly. Instead, the 
mere division of the functional relation (|26p at A = —E^ by its first conjugate 
partner leads to the pair of formulae 

2argDW±(-e-'^E,) - ^ p.,{v) = 7T[k + ^ ± ^] for k , (30) 
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which have the outer form of Bohr-Sommerfeld quantization conditions but 
are exact, albeit implicit (the left-hand sides invoke the determinants of the 
first conjugate spectra, equally unknown). 

However, we can write the totality of such conditions for the L conjugate 
problems (Sec. 2.1), together with the structure equations that express all the 
determinants from their spectra: namely, for i G Z/LZ and k,K 

logD[^+i]±(-e-'^Ejfl ) -logZ}[^-i]±(-e+i^Ef )j - {-lYip/3_i{v) (31) 
logDM±(A) = hm (j]log(4nA) + ilog(Ef + A) (32) 

{q>0} 

As a result, the exact spectra now appear as solutions of nonlinear, parity- 
split, selfconsistent or fixed-point equations of the form 

M+{S'+} = S*+, resp. M-{S-} = S- , (33) 

where denote the disjoint unions of the conjugate spectra f^'^^^. 

Numerical tests, done for < 6 and not too big coefficients v, suggest 
that the fixed-point equations (j!TT|) - (p?^ can be cast into contractive form. 
Assuming this conjecture, the problems (jH^ get solved simply by forward 
iterations, which induce convergence toward the exact spectra starting from 
trial spectra that need only be seniiclassically correct (to a certain order). 



4.3 Example: homogeneous potentials [144] — 

This is the simplest case (all conjugate spectra coincide, and A is the only 
variable), and specially so if A^ > 2 (see Sec. 4.1 for N = 2, jlH] for A^ = 1). 

For N > 2, the exact quantization conditions ()3H) reduce to a single pair 
(e.g., see fig.l for A^ = 4), 

2argD±(-e-'^E,)=7r[A; + i±^] for = JJJ;;;; (34) 
while the structure equations (jH^ simplify as 

argD±(-e-^'^E) = ^ arg(Efc -e-^'^E) (^ > 0), (35) 
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Figure 1: The exact quantization formulae (-I-, left; — , right) for the homo- 
geneous quartic potential g^: the curves 

{D^{e'''/^ E)}e>o (solid lines labeled by 
ii^- values) cross the (dashed) lines of prescribed argument in the complex D-plane. 



K , = (A'-2) / (A'+2) 




1 2 3 4 6 8 10 20 50 ~ 

N 



Figure 2: Numerically measured linear contraction factors k± at the fixed points 
(i.e., the exact it spectra), upon iterations of the exact quantization conditions 
(|34|) -(|35 |) for homogeneous potentials \q\^ ■ As horizontal scale we use the value 
taken by n± classically (when is replaced by as in ea. H21() ): Kci{N) = 



11 



and the asymptotic boundary condition is b^E^^ ~ A; + | for A; +00. 

The latter fixed-point equations have now been proven globally contrac- 
tive on suitable sequence spaces 0, by methods specific to the homogeneous 
case (see fig.2). (For the more general equations (jHT|) - (jH^ . contractivity 
remains a conjecture.) 

Those very same equations have also been found to yield the 

ground states of some exactly solvable models in 2D statistical mechanics 
[211 EH] and in conformal quantum field theories {Bethe Ansatz equations). 
A deep explanation for this coincidence of equations (and, by uniqueness, of 
their solutions) in so different dynamical settings is not yet available. This 
"ODE/IM correspondence" (as a dictionary between Schrodinger Ordinary 
Differential Equations and Integrable Models) currently applies to homoge- 
neous problems (even singular ones) plus the supersymmetric equations of 
Sec. 4.1 on the Schrodinger side, and only to the ground (non-excited) states 
in integrable models. There is now no integrable-model realization of the gen- 
eral inhomogeneous quantization conditions (j3T| - (j32|) : by extension, we still 
consider these as Bethe Ansatz equations, and the general ID Schrodinger 
equation as exactly solvable through them. 

5 Extensions: realized vs desirable 

5.1 Solving for the unknown functions i/; |47] — [^T| 

We keep focusing on the recessive solution ip\ of eq.((H) as an example. By 
letting now the reference half-line vary as [Q,+oo), we obtain parametric 
translates of the basic identities (j21|), essentially ipxiQ) = -Dq(A), ij'\{Q) = 
—Dq{\). The computation of V'a(Q) thus amounts to a parametric spectral 
problem, and this is now exactly solvable through parametric fixed-point 
equations like e.g., ^AqISq'} = £q~ for ip\{Q). (This has also under- 
gone successful numerical tests.) 

5.2 Toward quantum perturbation theory jHj, I52| 

By rescalings a la Symanzik, the perturbative regime for an anharmonic po- 
tential like V{q) = + gq^ amounts to the v —>■ +00 regime for the potential 
-|- vq"^. This limit is not manifest, and is actually singular, in our exact 
quantization formalism. We recently obtained explicit large-u behaviors for 
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such spectral determinants, basically as 



det^( 



dg2 



+ g^ + V' + A) det^( 



dg2 



+ + vq^' + A) 



(36) 



det^( 



^ + vq^i + A) 



det,^( 



^ + vq^i + A) 



(for V — * +00, with > M); i.e., the singular (divergent) large-v behavior 
is entirely concentrated in the classical parts of the determinants; these now 
express in terms of improper action integrals (cf. Sec. 3.1), 



and the latter can be estimated in the end [5^ . 
5.3 Generalizations? 

It is most desirable to push exact solvability beyond polynomial potentials. 
According to some exact results, rational potentials [21 E] and trigonometric 
polynomials jSSl UH] appear as realistic targets, as well as some differen- 
tial equations of order > 2 [H HH |22l 0] • The general Heun class is also a 
possible direction [HIj, which includes other important equations (Mathieu, 
Lame, . . . ). 

Of course, an even greater challenge is to tackle higher-dimensional Schro- 
dinger-type problems, in line with the original Balian-Bloch surmise itself 
[Bj. Exact asymptotic methods are still in the exploratory stage there, see 
for instance [111 EHl GDI E3 El • 

Therefore the subject is far from being closed! 
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